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Abstract. The renormalization of the third moment of the twist-2 flavour non-singlet Wil- 
son operator is given to two loops in the RI'/SMOM renormalization scheme. This involves 
renormalizing the operator itself and the two total derivative operators into which it mixes by 
inserting it into a quark 2-point function and evaluating the amplitudes at the symmetric sub- 
traction point. The corresponding two loop conversion functions are derived from which the 
three loop Landau gauge anomalous dimension is deduced. The full set of amplitudes for the 
two loop Green's function for each of the operators are given in both the MS and RI'/SMOM 
schemes. 
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1 Introduction. 



Recently a renormalization scheme has been introduced which aids the extraction of accurate 
measurements of matrix elements or Green's functions from the lattice for strongly interacting 
fields, [HEIE]. The scheme is named RI'/SMOM which stands for the modified regularization 
invariant (RI') scheme at the symmetric subtraction point where the infinities are removed in a 
fashion akin to the momentum subtraction scheme, [1]. The shorthand for the latter is MOM. 
Whilst this is a precis of the syntax it in effect precisely defines the method for performing the 
renormalization. The RI'/SMOM scheme is an extension of the original regularization invariant 
(RI) scheme and its modification (RI'), [5,i6j. These were developed earlier to aid the matching 
of lattice computations with the high energy behaviour determined in conventional perturbation 
theory. In essence those schemes are such that 2-point functions are rendered finite by choosing 
the renormalization constant so that at the subtraction point there are no 0{a) corrections 
where a = g'^ /{IGTr"^) and g is the coupling constant in QCD. By contrast the renormalization 
of 3-point functions such as vertices or operators inserted into 2-point functions is carried out 
using the standard MS scheme, [5l [6j. The technical difference between RI and RI' resides in 
the way the quark 2-point function is treated prior to applying the renormalization condition. 
More specifically this is how to project out the contribution associated with |^ where p is the 
external momentum of the 2-point function. The scheme which is more widely used of the two 
is RI' as it minimizes the number of derivatives required to be taken. This reduces the financial 
penalty in respect of defining the scheme when lattice regularization is used. By contrast using 
dimensional regularization in the continuum there is no such cost, only the limitation in the 
range of validity of the coupling constant, due to its smallness, which is not an issue on the 
lattice. There the nonperturbative structure can be fully explored. Further, there have been 
several continuum computations performed in perturbation theory in QCD, [71 [8]. For instance, 
the three and four loop renormalization of QCD in RI' was carried out in the Landau gauge in 
[7] and for an arbitrary linear covariant gauge at three loops in [8]. In [8] all 2-point functions 
were made finite in the same way as the quark 2-point function was in the original article, [5l[6]. 
Though it should be noted that in the Landau gauge the coupling constant and gauge parameter 
variables are the same in MS and RI'. 

We have given a resume of the RI' scheme partly to contrast with that of RI'/SMOM which 
is the main focus here, but also since it is retained as part of the latter scheme. This is because 
RI'/SMOM differs from RI' in the way 3-point and higher functions are treated. For instance, it 
was pointed out in \T] that RI' is sensitive to infrared effects. This stems from the fact that there 
is an exceptional momentum configuration for 3-point functions and the nullified momentum 
fiowing through an inserted operator leads to unwanted infrared singularities. Examples of 
where one has to be careful in the treatment are the quark current operators and the Wilson 
operators used in deep inelastic scattering. Therefore, there is a potential problem with making 
accurate measurements of matrix elements involving such operator insertions on the lattice which 
is the reason why RI'/SMOM was developed, [1]. In this scheme there is a momentum fiowing 
through each of the external legs and inserted operators of the Green's function in a way which 
is consistent with energy-momentum conservation. Therefore in the absence of an exceptional 
momentum configuration there ought to be no infrared problems. Indeed it was shown in [21 [3] 
that there was an improvement in the convergence of the conversion function associated with 
the renormalization at two loops when the RI'/SMOM scheme was used instead of RI' for the 
case of the scalar and tensor quark currents. More recently the second moment of the twist-2 
flavour non-singlet Wilson operator was studied in [9] at two loops. This built on the initial one 
loop calculation of [lOj and the two loop computations for the quark currents, [2l[3l[TT]. The 
Wilson operator was a more involved computation than the quark currents due to the mixing of 
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the operator with a total derivative operator. The latter is completely passive in the RI' scheme 
renormalization as the operator insertion is at zero momentum by definition. For RI'/SMOM 
it cannot be neglected and the mixing, which has been studied at three loops in MS in |12j . 
is crucial. Though comparing similar conversion functions for this operator moment suggested 
that, if anything, in the RI' case the perturbative convergence of the series was marginally 
quicker, [9]. Though it was noted in [9] that as the RI' scheme does not access the off-diagonal 
part of the mixing matrix, due to the form of the Green's function momentum configuration, 
it is in some sense not a full scheme for operators with mixing, in contrast to multiplicatively 
renormalizable operators such as the quark currents. 

Whilst we have concentrated on the renormalization, one ingredient which is crucial for lattice 
computations is the actual structure of the Green's function with the operator inserted in the 
quark 2-point function. This is used to measure nonperturbative matrix elements. However, in 
order to assist matching the lattice results accurately with the continuum, the provision of the 
same quantity in perturbation theory to as high a loop order as is computationally possible is 
important. Indeed the RI'/SMOM lends itself naturally to this problem of measuring forward 
matrix elements and hence generalized parton distribution functions. For instance, there has 
been significant lattice work in this general area represented by [l3l El [ISl [ISl [171 [IHl [IHl EQl 
[211 [221 [231 [Ml [251 lig [271 [281 [29]. As the second moment of the flavour non-singlet Wilson 
operator has been treated in [9] , it is the purpose of this article to extend that work to the case 
of the third moment building on the one loop result of [lOj. This is a larger computation as 
the extra covariant derivative means that there is mixing into two total derivative operators, 
[12]. However, as discussed there these operators are actually related to the total derivative of 
both second moment operators. In practice this means that the anomalous dimensions should 
be the same but only some of the amplitudes of the Green's functions will be related. This 
latter property is due to the Lorentz index imbalance between operators of different moments. 
In addition to focusing on the RI'/SMOM scheme structure of the Green's function we will also 
provide the full set of MS results at the symmetric subtraction point. Aside from being the 
reference renormalization scheme we record these because some lattice groups choose to devise 
their own RI'/SMOM type scheme and then convert their results to the MS scheme. Therefore, 
they would perform the high energy matching to perturbation theory in the continuum in MS. 
Whilst all the lattice computations are performed in the Landau gauge, our results will be in an 
arbitrary linear covariant gauge. This is because the presence of the gauge parameter is partly 
used as a checking procedure within the symbolic manipulation programmes we have used. In 
addition as was noted in [9l [TO] when one has two or more free Lorentz indices on the inserted 
operator, then there is not a unique way of defining RI'/SMOM. This is mainly to do with 
what combination of basis tensors one uses to write the Green's function itself in. Different 
basis choices could lead to different RI'/SMOM schemes. Indeed it might be possible to choose 
one in such a way that the associated conversion functions have a better convergence than the 
analogous RI' one. One final point concerning the third moment of the Wilson operator and that 
is that the vector current is present in the set of operators at this level through one of the total 
derivative operators, [12j . As this current is physical then its renormalization is determined by 
general properties of renormalization theory and the Slavnov- Taylor identities. Therefore, we 
will pay special attention to that to ensure there is consistency. 

The article is organized as follows. We devote the next section to the background to the 
set of operators we will renormalize and the structure of the Green's function we will evaluate 
to two loops in both schemes as well as a summary of the technical machinery we employed to 
perform the computation. The details of the two loop renormalization as well as the RI'/SMOM 
anomalous dimension mixing matrix are given in section three. The explicit forms of all the two 
loop conversion functions are given in the subsequent section. These are then used in section five 
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to determine the Landau gauge anomalous dimensions for the diagonal and some off-diagonal 
entries of the mixing matrix at three loops. The expressions for the amplitudes of the Green's 
function we compute are given in section six for both the MS and RI'/SMOM schemes with 
conclusions provided in section seven. An appendix gives the explicit form of the tensor basis 
used for the Green's function as well as the projection matrix used to extract the amplitudes. 



2 Background. 



We begin by recalling the formalism and properties of the operators which we will consider here. 
The notation we use is the same as in previous articles, [9l [lOl [TTl [12], and we will denote our 
three operators with the shorthand notation. 
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where all derivatives act to the right and S denotes both symmetrization and tracelessness, in 
d spacetime dimensions, in all Lorentz indices. In particular 
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As in previous articles the symbol VF3 at times will be referred to as the level and in those 
instances, which will be clear from the context, indicates the set of the three operators. For 
example, it appears here in the superscript. The specific choice of level VF3 operators is dictated 
by the need to have three independent operators. One could have excluded one of the total 
derivative operators in favour of one with the covariant derivative acting on the anti-quark field 
but we retain consistency with lower moment operators here, [9l [TOl 112) . In ()2.ip we have 
suppressed the flavour indices and emphasise that we are using flavour non-singlet operators. If 
the operators were flavour singlet then one would require gluonic and ghost operators to have 
a closed set under renormalization. In the flavour non-singlet case there is mixing but only 
within the set ()2.ip which is also closed under renormalization. However, the mixing is only 
relevant when the parent operator of W3 is inserted in a Green's function with a momentum 
flowing through the operator itself, [12]. With the way we have chosen the set of operators then 
if no momentum flows through the inserted operator there is no contribution to the Green's 
function from the total derivative operators. For the RI'/SMOM scheme renormalization, [1], 
where there are two external momentum, p and q, flowing in through the quark legs then there 
is a net momentum flow of [p + q) out of the operator, as illustrated in Figure 1 where the form 
of the Green's function, (^l;{p)0^^y„{—p — q)'tp{q)'^, we consider is given. Since we work at the 
symmetric subtraction point, [Il[2l[3], then 
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Figure 1: The Green's function (il){p)0^ {—p — q)i^{q) 



where /x is the common mass scale which in this case is equivalent to the mass scale introduced 
to ensure that the coupling constant is dimensionless in dimensional regularization which we use 
throughout. 

is accommodated by the mixing matrix of 
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The renormalization of the operators of (12.1 
renormalization constants defined by 
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where the subscripts label the operators in the order W3, dW^ and ddW^,. The upper triangular- 
ity of the matrix follows from the choice of our operators in (j2.ip and simplifies the computation 
when compared to the extraction of all nine elements for another choice. Associated with the 
renormalization constants are the anomalous dimensions which are encoded in a parallel matrix 
with the formal definition 
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Here a is the shorthand for the coupling constant via a = ^^/(IGvr^) and a is the canonical 
gauge parameter associated with a linear covariant gauge fixing with a = corresponding to the 
Landau gauge. Further, /3(a) is the /3-function and 7a(a, a) is the anomalous dimension of the 
gauge parameter where we retain the conventions of [8]. Whilst the anomalous dimensions for 
gauge invariant operators are independent of the gauge parameter in massless renormalization 
schemes, such as MS, we have retained it here since the RI'/SMOM scheme is a mass dependent 
scheme and therefore the anomalous dimensions will depend on the choice of gauge as will 
be evident later and hence is indicated in ()2.8p . The explicit values for the MS anomalous 
dimensions have already been determined in the MS scheme and we record them as they are 
required for constructing the RI'/SMOM three loop anomalous dimensions later. From [12] 
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where C(z) is the Riemann zeta function and the colour group Casimirs are defined as 
T^T" = Cf , /'^^'^/^^'^ = CaS"^ , Trfr^r'') = S^^Tf 
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where 1 < a < Na, Na is the dimension of the adjoint representation and Nf is the number of 
massless quarks. However, only the diagonal elements and the 23 elements are known to three 
loops. The off-diagonal elements for the first row were only determined to two loops. The reason 
for this is that to extract the mixing matrix renormalization constants, |12j . one has to consider 
various momentum routings of the operator inserted in the Green's function of Figure 1 as well 
as exploit the fact that there must be no terms such as {ln{p^ / fi'^))/e, where d = 4 — 2e and p 
is the sole momentum flowing through external legs. This latter criterion establishes relations 
between the off-diagonal counterterms. For the n = 2 moment operator one had sufficient linear 
equations in order to determine the three loop mixing matrix, [12j. With the increase in moment 
one requires a further linear relation in order to resolve the counterterms at three loops. This 
can only be achieved by a four loop renormalization for which was not considered in [12j. 

Throughout we will use the convention that the results will be expressed in the RI'/SMOM 
scheme, using RI'/SMOM variables, unless otherwise specified as shown in (I2.10p . As we will be 
mapping between the MS and RI'/SMOM schemes we recall that the variables are not necessarily 
the same in each scheme. However, the relations are known to three loops, [8], and are 
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which imphes that the Landau gauge is preserved under mappings between schemes. These 
relations derive from [8j where the full three loop renormalization of QCD in a linear covariant 
gauge in the RI' scheme was recorded. That scheme was defined by ensuring that all the 2-point 
function renormalizations were performed in such a way that there were no 0{a) corrections at 
the subtraction point. So included within this is the gluon 2-point function which determines 
the renormalization of a. The coupling constant renormalization was performed using 3-point 
function renormalization but those functions and higher are renormalized in the MS fashion 
whence the trivial mapping of (I2.12p . The RI'/SMOM scheme differs from RI' in that the 
3-point and higher functions are rendered finite by ensuring there are no 0(a) corrections at 
the symmetric subtraction point. We should note that in [21 13] the renormalization of a was 
assumed to have been carried out in the MS scheme. Whilst this means that Landau gauge 
expressions will be similar, if one was comparing the a dependent parts for the same quantities 
in the definition of [H [21 [3] and that used here, then they would not be in agreement for non-zero 
a. 

The Green's function we will compute with the level operators inserted is illustrated in 
Figure 1 where there are three Lorentz indices associated with the operator insertion which is 
denoted by the circle containing the cross. Therefore, at the symmetric subtraction point the 
Green's function has to be decomposed into a basis of Lorentz tensors which obey the same 
tracelessness and total symmetry as the original operator. The choice of basis tensors was given 
in [10] for the one loop computation and we retain the same basis here. More explicitly at the 
subtraction point we formally write 
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where T''(^k) fiua^P^l) basis tensors. For there are fourteen of these and as they are 

cumbersome we have provided the explicit forms in the Appendix. It should be noted that 
the tensor basis we use is defined at the symmetric subtraction point and is the same basis for 
each of the operators of level W3. When the squares of the external momenta of the Green's 
function are not all equal then the basis of tensors will be much larger. The objects which 
we will concentrate on computing are the amplitudes associated with each tensor for each W3 
operator insertion. These are T,'jj^^{p,q) and defined by the projection method, [TO] . 
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Here M\i is a dimension fourteen matrix of rational polynomials in d. It is computed by inverting 
the matrix Af^^ which is defined by 
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Due to the large number of basis tensors this matrix is equally cumbersome and is also recorded 
explicitly in Appendix A. 

In the choice of basis tensors we used the generalized 7-matrices F^^y'^"', [301 [3T1 [32 l [331 IM] . 



which are defined by 
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where the notation includes an overall factor 1/nl. These are totally antisymmetric and form 
the generalization of the 7-matrices in d-dimensional spacetime. Athough we are ultimately 
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interested in four dimensions they are a more appropriate object to use to build the tensor basis 
as the generahzed 7-space naturahy partitions into regions defined by the number of free Lorentz 
indices. The justification for this hes in the observation, \30 \ I3 H [32]. 



where /Mi - Mmi'i---fn jg ^j^g ^j^j^ matrix in this F-space. Indeed this is why (1A.2P has a F^-j^^ 
and ^^^-j^ subspace. Though we will always use 7^* as the object rather than the generalized 
object. Hence, when constructing the tensor basis one can be confident that a complete set 
has been used. For instance, the antisymmetry means that at most two external momenta 
can be contracted with F^^j"'^" for n > 2. So the Lorentz aspect of the basis is built from 
combinations of tensors from the set {r]^'^ , p^^ , ,T^^^"^"} . The choice of the fourteen given in 
Appendix A were also derived from the explicit one loop computation, [10]. In other words 
the one loop Feynman diagrams were expressed in terms of the basis Lorentz tensor Feynman 
integrals which were evaluated directly. The expressions for these were substituted back into the 
original computation and the Lorentz indices contracted producing the fourteen basis tensors, 
[lOj. Therefore, given that we use the basis here but do not construct all possible two loop basis 
Lorentz tensor master integrals. Instead we use the projection method. 

The machinery used to perform the computations are the symbolic manipulation language 
Form, [35] . The 3 one loop and 37 two loop Feynman diagrams are generated in electronic 
format using the Qgraf package, [36], and then Lorentz and colour group indices appended. 
As the Green's function contains two external momenta then packages such as Mincer, ^37J, are 
not applicable. Instead we first rewrite all the diagrams as scalar integrals where the numerator 
tensor structure is rewritten in terms of the denominator propagators in as far as this is possble. 
In this form we then used the Laporta algorithm, [38]. This creates a set of linear relations 
between all integrals with a certain specification for each appropriate momentum topology. The 
relations are established by integration by parts and Lorentz identities and the integrals of 
interest are rewritten in terms of a relatively small set of basic scalar master Feynman integrals. 
For the two loop Green's function these have already been established to the order in e we need 
in order to find the finite part. This is not an empty statement because in the huge set of linear 
equations spurious poles in e can emerge in the factor multiplying the master integrals. To handle 
such an enormous set of equations requires computer technology. Several packages are available 
but we used the Reduze programme, [39], based on GiNAC, |40j . which is written in C++. It 
turns out that for the 37 two loop diagrams there are only two basic momentum topologies. In 
other words the momentum routing of all the diagrams could be expressed in terms of one or other 
of these configurations. Also given that we had constructed the algorithm for the lower moment 
operators, W2, it was straightforward to lift the extra integral relations needed for W3 out of the 
earlier database created for W2- Therefore, we are reasonably confident that our procedure is 
consistent as the same routines were used for the quark current operator renormalization which 
agreed with the two loop results of [21 [3] for the scalar and tensor operators. 

3 Renormalization. 

The renormalization of the W3 operators has been outlined in [12] and we briefiy recall the 
procedure. First, we compute the Green's function with all the parameters being regarded as 
bare. This is the simplest way to proceed when dealing with automatic symbolic manipula- 
tion computer programmes. The renormalized variables are then introduced via the associated 
renormalization constant in the canonical way, [41j. This ensures that the counterterms are cor- 
rectly embedded within the Feynman diagrams automatically. For the operators themselves the 
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renormalized versions are also introduced by a similar rescaling involving their renormalization 
constant. However, unlike the fields and parameters this renormalization is not multiplicative 
due to the mixing, [12]. Therefore, when the renormalized operators are introduced we include 
the operators into which it mixes. Consequently poles in e will appear in various channels of 
the tensor basis but not all. They have to be removed by the criterion of the scheme in ques- 
tion. For MS this is achieved in the standard way. However, for the RI'/SMOM scheme it is 
more involved. If, for the moment we concentrate on a generic tensor channel which contains 
divergences, and denote that channel by 0, then the renormalization condition is, [1], 



The quark wave function renormalization is performed in the RI' scheme but as we are now 
dealing with a 3-point function then the operator renormalization constant obeys the same 
condition. The ethos for the regularization invariant schemes is that after renormalization there 
are no 0{a) corrections at the subtraction point which here is the symmetric point, [Ij. So 
both the quark 2-point function and the Green's function we are considering have no 0(a) 
corrections. For W3 we need to qualify ()3.ip by saying that the generic renormalization constant 
could represent a combination of counterterms which appear in the mixing matrix. Therefore to 
determine all the information to compute the anomalous dimension mixing matrix requires the 
solution of a set of linear equations. This was also a feature of the original MS renormalization of 
the operators in [12]. Whilst (13. Ih is the condition which determines the operator renormalization 
constant to avoid confusion it is important to note that the full Green's function is multiplied 
by the mixing matrix of renormalization constants. So after the renormalization constants have 
been fixed the finite parts of all the amplitudes are affected by the renormalization. 

Having outlined the general aspects of the renormalization we have to discuss the specific 
renormalization of ddW^ which requires special attention. This is because that operator is 
related to the divergence of the vector current. The vector current is a special operator in that 
it is physical. Therefore, not only does it not get renormalized to all orders in perturbation 
theory but the renormalization constant is the same in all schemes. See, for example, [12] for a 
summary of this. This is a consequence of gauge symmetry and effected by the Slavnov- Taylor 
identities. Therefore, our RI' scheme renormalization of ddW^, must be consistent with these 
general principles. So from the definition of ddW^ we have to project out that piece which 
corresponds to the total derivative of the divergence of the vector current. Contracting the 
Green's function containing the ddW^ operator with (p + q)^riua does not produce anything 
non-trivial due to the traceless condition. Instead we contract with {p + q)fj,{p + q)u{p + Q)a 
and this combination of vectors is chosen since that is the momentum fiow through the operator 
itself. This produces a non-trivial combination of the amplitudes. Ordinarily this would produce 
an expression where there are poles in e. Indeed the individual amplitudes have divergences. 
However, as the underlying operator is finite there are no divergences in e. Whilst this is 
consistent with earlier remarks one has also to ensure that the Slavnov- Taylor identity is actually 
satisfied and this relates to the finite part of this projection. Computing in the MS scheme, where 
the scheme dependence is in effect the quark wave function renormalization from (|3.ip . we find 
that the combination of amplitudes of the contraction proportional to |i is, |10] . 
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A similar expression for the piece involving ^ produces the same outcome. Clearly it is propor- 
tional to the finite part of the quark 2-point function after renormalization in the MS scheme. 
Indeed the equality of the ^ and ^ parts with the unit renormalization constant required for 
ddWs, due to the absence of poles in e, as well as the explicit values for the MS amplitudes 
means that our renormalization is consistent with the Slavnov- Taylor identities. Repeating the 
process for the RI'/SMOM scheme with the RI' scheme quark wave function renormalization 
constant produces a similar result which is 
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Clearly there are no 0{a) corrections which is in agreement with the RI' scheme quark 2-point 
function after renormalization. The fact that when this combination is computed that there 
were no poles in e means the renormalization of ddW^ also has a unit renormalization constant 
in this scheme consistent with the physicality of the operator itself. 

The procedure to determine the remaining renormalization constants of the mixing matrix 
is to identify a set of amplitudes or their combinations to which one can apply the condition 
(j3.ip . As noted in [9l [10] for the RI'/SMOM scheme there is not a unique way of doing this. 
One could use the same combinations for W3 and dW^ as ddW^. However, as there are five 
renormalization constants to determine one would require four other conditions. Instead we 
extend to two loops the RI'/SMOM scheme definition given in [iQj. Here the first three channels 
are rendered finite by the generic condition (j3.ip . In addition to the lack of uniqueness in defining 
the scheme because of the choice of amplitudes there is also the issue that our basis choice of 
tensors defining the amplitude channels is not unique. One could have defined a different set of 
basis tensors, though with the same symmetries as the operators themselves. Then one would 
have another definition of RI'/SMOM. The choice we make here is in some sense a minimal 
choice and one could readily make other more exotic ones. However, following this particular 
path we have determined the renormalization constants and have encoded them in the mixing 
matrix of anomalous dimensions. We have 
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+ 



1152V''(i) - 7687r^ 



24696 



Ca - 6105Cf 



TpNf 



1296 



7i^'(a,a) 



RI'/SMOM 



-Cpa 



+ 



(2112 - 864a - 288a^)4^'{^) + (1920^ + 576a - 1408)71^ 
+ 216a2 + 648a - 636841 Ca + 2187Cf 
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RI'/SMOM 



(a, a) 



RI'/SMOM 



(a, a) 



RI'/SMOM 



7^' (a, a) 



RI'/SMOM 



+ 



26280 + 5127r2 - 768tp'{l) TpNj 



3888 



+ 0(a^ 



-Cpa 
2 



(288a2 + 864q - 1056)^'(i) - (1920^ + 576a - 704)7r2 



1188a^ - 3564a - 10872 



Ca 



+ 



+ 



(2592 - 2784a)V''(3) + (1856a - 1728)7r2 
8244a - 9207] Cp 



384V''(i) - 2567r^ + 2952 TpNf 



3888 



+ 0{a' 



-Cpa 
3 



;i08a^ + 324a- 924)V''i 
- 81a^ - 243a + 16866 



(72a2 + 216a - 616)7r2 



Ca 



+ 



336V''(5) - 2247r^ - 5976 TpNf 



2016Cf 



486 



+ 0(a3 



-Cpa 



+ 



176tt^ - SW - 243a - 6651 



Ca 



+ 
+ 

0{a') 



144^'(i) _ 288Tp'{l)a - 288 + 648a - 967r2 + 1927r2a^ Cp 

Cpa^ 



647r^ + 2340 - 96V''(i) TpNf 



486 



+ 0{a') , 



(3.4) 



where ip{z) is the derivative of the logarithm of the Euler F-function. As a check on our proce- 
dures we first verified that the two loop MS mixing matrix of anomalous dimensions emerged and 
agreed with [12] . Whilst [12] was a three loop computation it used the Mincer algorithm, |37j . 
designed for massless 2-point function renormalization in dimensional regularization. The one 
loop terms are clearly gauge parameter independent but the two and higher loop corrections will 
depend on the gauge parameter a. This is because in mass dependent renormalization schemes, 
such as RI'/SMOM, the anomalous dimensions of gauge invariant operators are not gauge in- 
dependent in contrast to mass independent schemes such as MS. A final observation on the 
calculation is that the renormalization of SWs is the same as that for W2 which was considered 
in [9] . This is a non-trivial check since one has a different set of tensors for the basis due to the 
imbalance of Lorentz indices between the operators. That the same renormalization emerges is 
partly due to a similar renormalization scheme choice but also because the set of operators Wn 
are all connected via the tower of operators involving the associated total derivative operators. 



4 Conversion functions. 

Having constructed the RI'/SMOM mixing matrix of anomalous dimensions at two loops we can 
construct the associated conversion functions. These allow one to translate between renormal- 
ization schemes. For us the two schemes will be RI'/SMOM and MS and general background 
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can be found in, for example, jl2]. However, as the renormalization of the operators is via a 
mixing matrix one has to extend the theory to ahow for this. So the natural extension is to a 
matrix of conversion functions which is formally defined by 



yW3 

fci.MS 



(4.1) 



As the parameters a and a are tied to a renormalization scheme we have to make a choice and 
note that in the argument of the conversion functions the parameters will be in the MS scheme. 
So in this definition those parameters in the RI'/SMOM renormalization constants have to be 
converted to the corresponding MS parameters via (|2.12p and (j2.13p . Otherwise aside from being 
inconsistent one would have a non-finite expression due to poles in e. For practical purposes, 
the definition translates into 

^W3( \ _ ii,Rl7SM0M „Wi( N _ 12,Rl7SMOM ll,Rl7SMOM^i2,MS 

^ii la'"J - -^wTi ' la'«J - yW3_ yWi vH-3_ ' 

M,MS 22,MS ll,MS 22,MS 

yW3 yWz yW-j, yW'i 

^W:i( ^ _ 13,Rl7SMOM ^ll,Rl7SMOM^i2,MS 23,MS 

^13 — „W3 7W3 yW3 yWs 

33,MS 11,MS 22,MS 33,MS 

yW3 yW3 yWs yW3 

ll,Rl7SMOM^13,MS 12,Rl7SMOM 23,MS 



yW3 vW3 yW3 vW3 ' 

ll.MS 33,MS 22,MS 33,MS 

yW3 yW3 yW3__ 

^W3i N _ 23,Rl7SMOM ^22,Rl7SMOM^23,MS „x 

(^23 - ^ty3_ Z^3_^H-3_ ^^■'^^ 

33,MS 22,MS 33,MS 



for individual entries of the conversion function matrix where the first term corresponds to the 
diagonal and there is no sum over i and i = 1, 2 or 3. Equipped with these the explicit values 
of the matrix elements to two loops are 

C{[3(o,q) = 1+ [(192a -216)^'(i) + (144 - 128a)^2_ 144^^ 2763] '^^'^ 



324 



+ 



3^ 

(294912a^ - 663552a + 248832) (^-'(1))^ 
+ (884736a - 3932160^ - 331776)V''(^)7r^ 

+ (13360896a - 857088a2 - 16819488)V''(|) - (5184 + 82944a) ?/'"'(i; 
+ (40310784a - 94058496)s2(f ) + (188116992 - 80621568a)s2 (f ) 
+ (156764160 - 67184640a)s3(f ) + (53747712a - 125411328)s3(|) 
+ (131072a2 - 73728a + 124416)7r^ 
+ (571392a2 - 8907264a + 11212992)7r^ 

+ 6168960^ - 3348864a + 20576997 + (1492992a - 1679616)S 

ln2(3)7r 



(3359232a + 3919104)C(3) + (279936a - 653184)- 



^/3 

Cf 



(7838208 - 3359232a) i^^^^ + (701568 - 300672a) ^ 

V 3 \/3 

62208('(/''(|))^ - 82944V'' (i)7r2 

+ (580608a2 - 1111968a + 974592)V''(i) 

+ (145152 - 46656a)V''"(|) + (65505024 - 25194240a)s2(|) 
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+ (50388480a - 131010048)s2(f ) + (41990400a - 109175040)s3(f ) 
+ (87340032 - 33592320a)s3(|) + (124416a - 359424)7r^ 

- (387072a2 - 741312a + 649728)7r^ 

- 505440a2 - 1714608a + 79566624 
(746496a - 1213056)S + (1679616a - 24354432)C(3) 



(454896 - 174960a) i^^-^ + (2099520a - 5458752) i^^^ 



Ca 



+ 

+ O(a^) 



- (187920a -488592)^ 
v3 

2052864V''(|) - 13685767r^ - 32160888] TpNf 



839808 



(4.3) 



(48 - A%a)^'{\) + (32a - 32)7r^ + 36a - 927 



+ 



324 



1290240^ - 396288) (V/(^))^ 



368640a + 528384)V'(i)vr 



(276480a 
+ (172032a2 

+ (4216320^ - 2847744a - 2817504)^'(i) - (155520 - 20736a)V'" 

- 73903104s2(|) + 147806208.S2(§) + 123171840s3(|) - 98537472s3(|) 

- (57344a^ - 67584a - 238592)7r'' 

- (281088a^ - 1898496a - 1878336)7r2 - 383616a2 + 440640a 

- 17134821 - (373248a - 373248)S + 13996800C(3) 



- 513216 



+ 



+ 



ln^(3)7r 



+ 6158592 



ln(3)7r 



IT 



+ 551232^ 
V3 



Cf 



124416(V''(i))^ 



165888V'' (5)7r2 



- (207360a^ + 147744a - 2208384)V''( 
+ (18144 + 15552aX'(i) + (15116544 + 5038848a)s2(|) 

- (10077696a + 30233088)s2(f ) - (8398080a + 25194240)^3 (f) 
+ (20155392 + 6718464a)s3(f ) + (6912 - 41472a)7r'' 

+ (138240a2 + 98496a - 1472256)7r2 + 225504a2 + 594864a 

- 23258016 + (93312 - 186624a)S + (839808 - 839808a)C(3) 

+ (34992a + 104976) ^^^^^^ - (1259712 + 419904a)-^''^^^'' 



V3 



V3 



vr 



- (112752 + 37584a) —= 
v3 



Ca 



3732487r2-559872V''(i) + 11611512] r^AT^]^^ + 0{a') , (4.4) 



(48a - 24)V''(|) + (16 - 32a)7r2 - 198a - 9 



(129024a^ - 193536a + 322560) (V-d)) 



Cpa 
324 



(172032a^ - 258048a + 430080)^/^ (g)7r 
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- (743040a2 - 1156032q + 9328608)V''(5) - (134784 + 20736aX'(i) 

- 80621568s2(f ) + 161243136s2(f ) + 134369280s3(|) - 107495424s3(f ) 
+ (57344a^ - 30720a + 502784)7r^ 

+ (4953600^ - 770688a + 6219072)7r2 + 974592a2 + 629856a 

- 2984931 + (373248a - 186624)S + 11757312C(3) 

ln2(3)7r ln(3)7r vr^ ' 

- 559872 y- + 6718464 y + 601344^ 

v3 v3 v3 



+ 



+ 



Cf 



165888V'' (|)7r^ - 124416(V''(|))^ 

+ (129600a2 + 1314144a + 1428192)^-' (5) 

+ (20736 + 15552a)V'"'(5) + (12597120 + 10077696a)s2(f ) 

- (20155392a + 25194240)s2(f ) - (16796160a + 20995200)s3(|) 
+ (16796190 + 13436928a)s3(f ) - (110592 + 41472a)7r'' 

- (86400a^ + 876096a + 952128)7r^ - 5171040^ - 1662120a 

- 1840968 - (93312 - 186624a)E + (629856 - 1679616a)C(3) 

+ (69984a + 87480) i^^-^ - (1049760 + 839808a) ^^^^^ 



7r3 

(93960 + 75168a) ^ 
v3 



Ca 

311040V''(i) - 2073607r2 + 390744] TpNf] + 0{a^) , (4.5) 



+ 



(36a - 42)V''(i) + (28 - 24a)7r2 - 27a + 4591 ^ 

J 81 

(5184a2 - 12096a - 4608)(V''(|))^ + (16128a - 69120^ + GIU)^)' {^)7r^ 

+ (328536a - 136080^ - 613656)^^' (3) - (5022 + 1944a)V'"'(^) 
+ (1259712a - 5248800)s2(f ) + (10497600 - 2519424a)s2(f ) 
+ (8748000 - 2099520a)s3(f ) + (1679616a - 6998400)^3 (f) 
+ (2304a^ - 192a + 11344)7r^ + (9072a^ - 219024a + 409104)7r^ 
+ 7290a2 - 90882a + 107568 + (34992a - 40824) S 

+ (104976a + 559872)C(3) + (8748a - 36450)—^ 

v3 

+ (437400 - 104976a) + (39150 - 9396a) ^ Cf 

v3 v 3 

+ [5832(V''(i))^ - 7776V''(|)7r2 + (11664a2 - 39366a + 99468)V''(|) 
+ (5103 - 972a)V'"'(|) + (2519424 - 629856a)s2(f ) 
+ (1259712a - 5038848)^2 (f) + (1049760a - 4199040)s3(f ) 
+ (3359232 - 839808a)s3(f ) + (2592a - 11016)7r^ 
- (7776a^ - 26244a + 66312)7r2 - 8748a^ - 34992a + 1759644 
+ (17496a - 34992)S + (26244a - 734832)C(3) 

+ (17496 - 4374a) i^^^^ + (52488a - 209952) i^^^ 
V 3 V 3 
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TT 



(4698a - 18792)^ 



Ca 



+ 



46656'i/''(3) - 311047r2 - 642168 TpNf 



26244 



+ 0{a^ 



(4.6) 



"-"23 



(a, a) 



24V''(i) - 167r2 - 54a - 297 



+ 



(13824a + 35136)(V''(3))^ - (18432a + 46848)V''(3)vr2 
- (342144a + 19440a2 - 824904)^'(i) + 18792V^'"(i) 



(10917504 - 1679616a)s2(f ) + (3359232 - 21835008a)s2 (f ) 
(2799360 - 18195840a)s3(f ) + (14556672 - 2239488a)s3(f ) 



(6144a - 34496)7r^ 



+ 



+ (12960a^ + 228096a 
- 188892 + 46656S 



37908a^ + 76788a 

(139968a + 1609632)C(3) + (75816 - 11664a) 



549936)71-2 
ln2(3)^ 



+ (139968a 



909792) i^^^ + (12528a 
v3 



81432)^ 



Cf 



+ 



31104V''(i)7r2 - 23328(^'(i))2 + (107892a - 5832a2 - 122472)V''(i; 



+ (972a - 9720)^- (i) + (944784a - 4513968)s2(|) 

+ (9027936 - 1889568a)s2(|) + (7523280 - 1574640a)s3(f ) 

+ (1259712a - 6018624)s3(f ) + (15552 - 2592a)7r^ 



+ 



(3888a2 - 71928a + 81648)7r2 - 21870a2 - 113724a - 2488482 
34992S + (1277208 - 52488a)C(3) 



+ (6561a - 31347) 



(33669 - 7047a) 



ln2(3)7r 



+ (376164 - 78732a) 



vr 



Ca 



+ 



V3 

2937671^ - 44064V''(i) + 946080 TpNf 



_CF(f_ 

104976 



ln(3)7r 

~7r 



(4.7) 



and 



C^'{a,a) = 1 + 0{a^) . (4.8) 

The final expression merely reflects the fact that ddW^ is related to the vector current which 
is a physical operator. The values for row 2 and 3 correspond to the values of the conversion 
function matrix for W2 given in [9j3. This is consistent with the Wilson operator renormalization 
being part of a tower of operators. These expressions include the usual transcendental type of 
numbers such as powers of vr and the Riemann zeta function of odd argument as well as rationals. 
However, given the nature of the subtraction point other classes of basic numbers arise. These 
include derivatives of the Euler -(/^-function and the natural logarithm as well as the polylogarithm 
functions through the specific function 



(4.9) 



*The expression for (a, a) corrects a typographical error in the presentation of the two loop term involving 
Cp of the corresponding equation for C|^^ (a, a) recorded in but which was correct in the attached data file of 
that article. 
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In addition a combination of various harmonic polylogarithms also appears which is 

s = nfi + -Hfg) (4.10) 

using the notation of [2]. These particular basic numbers derive from the explicit expressions for 
the scalar master integrals which the Laporta algorithm produces as a result of the reduction of 
integrals. They have been evaluated by various authors in \A?>\ HS} 



For the conversion to all the numerical forms we present here we have used the following 
numerical values for the various polygamma and polylogarithm functions which are 

C(3) = 1.20205690 , S = 6.34517334 , = 10.09559713 , 

ijj'" = 488.1838167 , S2 = 0.32225882 , S2 = 0.22459602 , 

S3 = 0.32948320 , S3 = 0.19259341 . (4.11) 

Expressing the conversion matrix in numerical form for the colour group SU (3) produces 

C]^3(a,a) = 1 + [2.1853715a + 8.2451607]a 

+ [9.95942160^ + 38.4241424a + 224.4265963 - 19.8007988iV/]a2 + O(a^) , 
C]^^(a,a) = - [0.5463429a + 3.1203238]a 

- [3.4486064a2 + 15.2459694a + 91.8037478 - 7.6549878iV/]a^ + O(a^) , 
CY^^{a,a) = - [0.1203238a + 0.3842826]a 

- [0.3835601a^ + 1.2332093a + 9.0543723 - 1.1782990iV/]o^ + O(a^) , 
CY^\a,a) = 1 + [1.6390287a + 5.1248369]a 

+ [6.5108151a2+ 23.1781730a + 132.6228486 - 12.1458110iV/]o2 + O(o^) , 
C^^^{a,a) = - [0.4444444a + 1.7499534]a 

- [2.1301455a2+ 7.2772405a + 49.1967063 - 5.0243682iV/]a2 + O(a^) , 
cZ\a,a) = 1 + 0{a^) . (4.12) 

Clearly there is a large correction at two loops for itself compared with dW^ or W2- As one of 
the motivations for developing the RI'/SMOM scheme was the hope that the convergence of the 
conversion function would improve we recall the numerical values for the conversion functions 
from RI' to MS. From [47] we have the numerical values in the Landau gauge 

CYi'ia,^) = 1 + 7.9259259a + [215.0853593 - 18.9809671iV/]a2 + 0{a^) , 
C^a^Q^Q) = 1 + 4.5925926a + [119.8268158 - 10.9794239iV/]a2 + O(a^) , 
CY^\a,Q) = 1 + 0{a') . (4.13) 

Given the fact the top entry is the key one it appears that the RI'/SMOM scheme has a slightly 
larger two loop correction compared to the RI' conversion function. A similar observation was 
made for W2, \W\ - However, as noted in that article it is not entirely clear whether one can truly 
compare these conversion functions. This is because the nature of the RI' scheme is such that it 
cannot access the off-diagonal elements of the mixing matrix. This stems from the momentum 
configuration of the underlying Green's function. Indeed in this respect it is not entirely clear 
whether one can regard RI' as a full renormalization scheme for the Wilson operators as a result 
of the mixing with total derivative operators. For operators where there is no mixing such as 
the quark currents then a comparison of conversion functions would appear more appropriate, 
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5 Three loop anomalous dimensions. 



One aspect of the conversion functions is that they can be used to determine the anomalous 
dimensions given knowledge of the anomalous dimensions in one scheme. Although wc arc 
working with mixing matrices here it is straightforward to extend the formalism to show that 
in our case 



7, 



7SM0M ("RI'' "RI') - ^^ik'' ("MS' "ms) '^Zm ('^MS' "ms) ('^/T' ("mS' "ms)) 

(^S''("MS'"Ms))"' (5-1) 



^d/j,^^* ("ms'"ms. 



where we now have to explicitly label the scheme the variables correspond to by the subscript. 
More explicitly for the three diagonal elements of the anomalous dimension mixing matrix we 
have 



7. 



ii,m'/SMOM 



(^RI'' «RI') = T^is ("Ms) - ("Ms) 1^ ("MS' «MS. 





"ms7^^ ("MS' «Ms) ^i:^ 1" ("MS' "Ms) (5.2) 

^"ms 



where again there is no sum over i and ? = 1, 2 or 3. The off-diagonal elements are more involved 
since 



7IIV/SM0m("RI''"RI') = 7,'^iTs("Ms)C'n'("MS'" 



*ms) 



P ("ms) Qa—^^^ ("ms'"ms 



*MS 



) ^ry ("mS'"ms) 

""ms 



MS'-MSy 0^. 



+ 7, 



llMS ("ms) ^is'' ("MS' "ms) 

("Ms) ^12^ ("ms' "Ms, 



22,MS 



+ 



+ 



Ci2^ ("ms'"ms) ^ ("ms) l^C'n'^ ("ms'"ms 

^12^ ("MS' "Ms) "Ms7^^ ("MS' "Ms) 



da 



— In C{[3 (a^, a^) [c^^ (a^, a^^)] , (5.3) 
•MS J 



^13,RI7SM0M l"RI''"RI'i 



^iThS ("ms) («'MS' "Ms) 



^ ^23Ms ("ms) ("MS' "ms) 



llJ5S ("ms) ^^"^ ("MS' "ms) 

'"Ms) 



+ ^sTk ("Ms) l"MS 

- /5 («Ms) ^;;^^i3^ ("MS' "Ms) 
^"ms 



— a 



Ms7q 



'^^ ("ms'"ms) ^13^ ("ms'"ms) 

""ms 
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+ '^n' («MS ' "Ms) P («Ms) ^9;— 1^ ^n' («MS ' 
+ ("mS' "Ms) "Ms">'a^^ ("MS' "Ms 



MS; 



d 



X 



da- 



MS 



In CYi (a^S' "ms) .^33^ ("^MS' "MS, 



+ 



and 



1. 



W3 



23,RI7SM0M ("RI''"RI') 



7. 



W3 



11,MS V^MS 



"'MS'"mS 



'^12 Ms ("Ms) '^11^ ("MS' "Ms 



"ms'"ms 



+ /5(«Ms) /5|='^i^'(«Ms'" 



MS 



+ "MS^a^^ («MS'"Ms) f)r, ^12^ ("MS'"MS 

•^"ms 



("^MS' "ms) ^ ("ms) ("^ms' "ms) 

/ \ Ms / 

*-"i2 v"-ms'"msJ "ms^o V'^ms'"ms, 
= l-<^(«MS'"Ms)1^2T^( 



_d 

da- 



"ms'"ms 



("ms'"ms) ^sa"* (■ 



°'MS'"mS 



(5.4) 



^23 Ms (°Ms) ^22'' («ms'"ms) 



MS'"mS 



/3 («Ms) Q^Cf,^ (- 

""ms 

Ms / \ ^ f^Ws ( \ 

"MS^a l°'MS'"Msj f)^ ^23 l°^MS'"Msj 

'^"ms 

^22,Ms V^'MsJ ^"i^ v"ms' "ms, 
^3^Ms (""ms) ( 



'^MS'"mS 



'^23' ("^MS' "Ms) ^ ("Ms) 
^23'^ ("-Ms' "Ms) "Ms'^'o^^ ( 



MS 

"Ms'"Ms 



-,1^3 

--22 Io-ms'"ms 



d 



da- 



US, 



InC"^^ (a^g,a^g) ^iS''' (aMS'"Ms). • (^•^) 



Using these we have first verified that the two loop arbitrary gauge anomalous dimensions, (I3.4p . 
are reproduced. However, equipped with the explicit values of the two loop conversion functions 
as well as various three loop MS anomalous dimensions we have the Landau gauge expressions 



7iT^(a,0) 



RI'/SMOM 



25 



+ 



70471^ - 1056V''(|) + 22068 



Ca - 2035Cf 



18 



+ 



+ 



2567r2 - 8232] TpAT^] + 0{a^) 



[[25344(^'(i))2 - 33792V''(|)7r2 + 4454V(i) 

+ 59136V'"'(5) + 26687232s2(f ) - 53374464s2(f ) 

- 44478720s3(f ) + 35582976s3(f ) - 1464327r^ - 296967r2 



494208S - 9066816C(3) + 42261846 + 185328- 



V3 



+ 



+ 



+ 



+ 



2223936 



ln(3)7r 



- 199056^ 
V3 



25344(V''(i))^ - 33792^/''(i)7r2 



3^,, 4907232V''(i; 
3 ; ^^^^^.^^.-^y-^j + 76640256S2 (f ) 
+ 63866880s3(f ) - 51093504s3(|) + 168967r*^ 
+ 32714887r2 - 684288E - 969408C(3) - 4991286 

- 266112^^ + 3193344^^^ + 285824^' 

V3 V3 



2112V'"'(i) - 38320128s2l 



^/3 



CaCf 



12288V/ (I 
- 21504V' 



- 9216(V''(i))^ + 899328V''(i) 

- 9704448s2(f ) + 19408896s2(f ) 
+ 16174080s3(f ) - 12939264s3(§) + 532487r^ 

- 5995527r2 + 179712S + 497664C(3) - 29036688 



,/// / 1 ■ 

3- 



- 67392 



ln2(3)7r 



+ 808704 



ln(3)7r 



+ 72384— 
V3 



V3 V3 
12288V''(5)7r^ - 9216(V''(i))^ + 1908864V''(i) 
+ 768V''"(|) + 13934592s2(|) - 27869184^2 (f) 

- 23224320s3(f ) + 18579456^3 (|) - 61447r^ 

- 12725767r^ + 248832S + 2529792C(3) - 2559504 

.3 



+ 96768^^ - 1161216^^^ - 
2027527r^ - 304128V'' (5) + 4517952 



103936^ 
T^Nf 



CFTpNf 



+ [1710720C(3) - 733515] C]p 



Cpa^ 
46656 



(5.6) 



+ 



+ 



3087r^ - 462V'' (i) + 84331 Ca 



+ 



[l68V''( 



1127r2 - 2988 



TpNf 



1008Cf 



243 



+ 0{a^ 



[[64152(V''(i))^ - 85536V''(|)7r2 + 862812V''(i) 

+ 56133V'"'(|) + 27713664s2(f ) - 55427328s2(f ) 

- 46189440s3(|) + 36951552^3 (f) - 1211767r^ - 5752087r^ 

- 384912S - 7243344C(3) + 25273296 + 192456 " y 

v3 
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- 2309472 



ln(3)7r 



- 206712—= 

\/3 



cl 



+ 



■ 89496(V/(5))^ - 5901984V''(i; 



119328V/ (i)7r^ 
- 55242<'(i) - 57736800s2(f ) + 115473600s2(|; 
+ 96228000s3(|) - 76982400s3(f ) + 1075367r^ 
+ 39346567r2 - 449064S + 3639168C(3) - 4833270 

.3 



400950i^^ + 4811400^^^ + 430650^ 



V3 V3 
+ [31104(V''(|))^ - 23328V'' ( I )7r2 + 251424V''(i) 

- 20412V''" (i) - 10077696s2(f ) + 20155392s2(f ) 
+ 16796160s3(f ) - 13436928s3(f ) + 440647r^ 

- 1676167r^ + 139968S + 1259712C(3) - 16603056 



CaCf 



69984 



ln2(3)7r 



- 839808 



ln(3)7r 



+ 75168^ 

v's 



VS ' """""" 

+ [32544(V''(|))^ - 43392V''(|)7r2 + 2227824V''(|) 

+ 20088V'"'(^) + 20995200s2(f ) - 41990400s2(f ) 

- 34992000s3(f ) + 27993600s3(f ) - 391047r^ 

- 14852167r^ + 163296S - 559872C(3) - 742608 

'V3 



CATpNf 



559872C(3) 

+ 145800i^^Sl^ - 1749600^ - 156600-""^ 



+ 



V3 ^ 
1244167r2 - 186624V'' (i) + 2423520 



Tj.Nf 



+ [1679616C(3) - 90720] C| 



39366 



+ 0{a' 



(5.7) 



4 



264V''(i) - 1767r2 - 665lJ Ca + [l44V''(i) - 967r2 - 288j Cp 
+ [647r2 - 96V''(i) + 2340] TpNf] + 0{a^) 

342144V''(i)7r2 - 256608(V''(^))^ - 1082808V''(|) 

- 106920V''" (I ) - 49653648S2 (| ) + 99307296s2 ( f ) 

+ 82756080s3(f ) - 66204864s3(f ) + 1710727r'' + 7218727r^ 

+ 384912S + 12369672C(3) - 37423134 - 344817^"^ ^^^^ 



Cl 



+ 4137804^^^^^^ + 370359^ 

+ [477504(V''(i))^ - 636672V'' (^TT^ + 7978176V''(3) 

+ 204768V''" (I) + 117993024s2(f ) - 235986048s2(f ) 

- 196655040s3(|) + 157324032s3(f ) - 3338247r'^ 

- 531878477^ + 653184E - 11897280C(3) + 367416 
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In^fSlvr ln(3)7r 
+ 819396 y 9832752- ^ ' 



^/3 



+ 



880092^ 
^/3 



CaCf 



+ 



93312(V''(i))2 - 124416?/;'(i)7r2 - 150336?/^'( 

+ 38880<'(i) + 18055872s2(f ) - 36111744s2(f ) 

- 30093120s3(f ) + 24074496s3(f ) - 622087r^ 

+ 10022471^ - 139968S - 1749600C(3) + 24312312 

+ 125388^^^^ - 1504656^^^ - 134676^1 C^T^iV/ 
\/3 \/3 \/3 

208896V'' (^vr^ - 156672(V''(i))^ - 3297024V''(i) 

- 75168<'(i) - 43670016s2(f ) + 87340032s2(f ) 
+ 72783360s3(f ) - 58226688s3(f ) + 1308167r^ 

+ 21980167r2 - 186624S + 3079296C(3) + 4039632 

^3 



In^ 3 W In 3 W vr 

303264 + 3639168 V + 325728—= 

^/3 ^/3 



CfTfNj 



+ 



T^Nf 



176256V''(^) - 1175047r2 - 3494016 

138240V''(i)vr2 - lO368O('0'(i))^ + 1537056V''(i) 

- 34992<'(i) - 1679616s2(f ) + 3359232s2(|) 
+ 2799360s3(|) - 2239488s3(f ) + 472327r'^ 

- 102470471^ + 139968S - 1399680C(3) + 729648 



11664 



ln^(3)7r 

~7r 



+ 139968 



ln(3)7r 

~7f 



+ 12528—= 
\/3 



cl 



Cpa^ 
157464 

(5. 



and 



73'rn«,0) 



RI'/SMOM 



(5.9) 



Again the final expression is the same as that for the vector operator. We note that we have not 
given the complete matrix at three loops. This is because the complete three loop MS matrix 
is not known. The 12 and 13 elements were not determined in [15] as there was not sufficient 
information to disentangle the relation between the two counterterms. As was noted in [12] only 
a four loop computation could resolve this. However, the diagonal elements as well as the 23 
element are now available in the RI'/SMOM scheme for comparison with the MS expressions. 
In numerical form these anomalous dimensions are 



7iT^(«,0) 
72'^^KO) 



RI'/SMOM 



RI'/SMOM 



RI'/SMOM 



5.5555556a + [161.5815415 - 10.6202306iV/] 



+ 



6275.4956875 - 993. 60723 17iV/ + 24.6390623iV/ + 0{a 



3.5555556a + [104.7024244 - 6.5770518iVf] a^ 



+ 



4010.9803829 - 624.8817671iV/ + 14.9653337iV/ 



1.7777778a - [46.3028612 - 2.74688257V/] a^ 
1692.0143513 - 265.3339715iV/ + 6.0846171iV/j + O(a^) , 
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RI'/SMOM 

for SU{2>). For comparison the corresponding MS values are 



= 5.5555556a + [70.8847737 - 5.1234568iVfl a 

MS ^ 



MS 



+ 



'/J 

1244.9136024 - 199.6373883iV/ - 1.7620027iVf2 



3.55555560 + [48.3292181 - 3.1604938/Vy] 



+ 



859.4478372 - 133.4381617A^/ - 1.2290809iV/ 



_ = - 1.7777778a- [24.1646091 - 1.5802469iV/] 



429.7239186 - 66.7190809iV^ - 0.6145405iV/ 



MS ^ 



(5.10) 



a" + 0{a^) , 



(5.11) 



which illustrates that the higher loop corrections are numerically smaller in the MS scheme. 



6 Amplitudes. 

We devote this section to the main results which are the explicit forms of the two loop amplitudes 
in both the MS and RI'/SMOM schemes. Given the large number of amplitudes which have to 
be recorded we have chosen to do this in numerical form for an arbitrary colour group. In [9, llj 
we also recorded the exact two loop expressions in the form of tables using the notation 



0\{22)J22) 



^.^^)) C^CAa' + (EC)'r«^''^) + O(a^) • (6-1) 



Here On"* correspond to a basis set of numbers which naturally arise in the computation. These 
include the transcendental type of numbers which have appeared in the conversion functions 
and anomalous dimensions earlier. In addition the set an^ also included the gauge parameter 
dependence. The labelling is chosen in such a way that I indicates the loop order and in addition 
at two loops a second digit is appended to reference the colour group Casimir associated with 
that set. The other quantities, c^j)'n \ the actual rational coefficients of interest. The 
main reason for recalling this notation here is that attached to this article is an electronic file 
containing all the coefficients for all the amplitudes in both schemes in their exact forms in 
terms of the basis of numbers On^ . Therefore this is intended to allow for comparison with other 
results such as the relation of dW^ and ddW^ with W2 and dW2 respectively. 

Therefore we now carry out the mundane task of recording the results. First, in the MS 
scheme we hav^ 

T.^Up,q) = - [0.1280942 + 0.0401079a] a 



MS 



3.5739633 + 0.5909329a + 0.1445330a^ - 0.3927663iV, 



^Attached to this article is an electronic file containing the explicit analytic forms for all the amplitudes 
presented in this section. 
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- [0.6481481 + 0.1311651a] a 

- 11.9214676 + 0.6190048Q; + 0.4924713a2 - 1.6685976Ar/] + 0{a^) , 

- 0.3333333 + [1.5801848 + 0.0617905a] o 

+ [23.3909330 - 0.0454454a + 0.3587242a2 - 2.8780596iV/] + 0{a^) , 
[0.3425460 + 0.0484212a] a 

+ [4.5342388 + 1.1292147a + 0.17032770^ - 0.2361911^/] + O(a^) , 
[0.4238476 + 0.0751598a] a 



+ 



4.5935423 + 1.2759751a + 0.26196550^ - 0.3660535A(^] + O(a^) , 



[0.5842793 + 0.1503196a] a 



+ 



4.8574744 + 1.7081275a + 0.52372040^ - 0.684:1872Nf] c? + O(a^) , 



[1.0677459 + 0.6648679a] a 



+ 



23.3776082 + 5.5517350a + 2.41085910^ - 2.1261604A(^1 c? + O(a^) , 
[0.2222222 + 0.1553757a] a 

+ [7.1027076 + 2.4806545a + 0.54272100^ - 0.1633838^^1 o? + 0{a^) , 
[0.2829270 + 0.2088529a] a 

+ [9.2546037 + 3.0300483a + 0.73834240^ - 0.2645584Ary] + O(a^) , 
[0.4773248 + 0.3374900a] a 



+ 



15.4402123 + 3.8373519a + 1.2238427a^ - 0.7240161iV/] + O(a^) , 



[0.5615110 + 0.1503196a] a 



+ 



7.6079600 + 1.7765023a + 0.54395530^ - 0.6061888A(^] + O(a^) , 

- 0.0607048a 

- 1.2759196 + 0.2443179a + 0.00505870^ - 0.0499049iV/l + 0{a^) , 

- 0.1214095a 

- 2.5537039 + 0.3188150a + 0.01011750^ - 0.1270706Ar/] + 0{a^) , 

- 0.4588995a 

- 8.6460262 - 0.5239987a + 0.0382416o^ - 0.7578066iV/] 

+ 0(a3) , (6.2) 

- [0.5833178 + 0.1481415a] a 

- 12.4616428 + 0.8098939o + 0.55362850^ - 1.6747894iV/] + O(a^) , 

- 0.1666667 + [0.2708217 - 0.0971989a] a 
1.0208053 - 0.5371502a - 0.3019999a^ + 0.0393764iVy] + O(a^) , 



+ 
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MS 



MS 



MS 



MS 



MS 



MS 



MS 



0.3333333 + [1.1249617 - 0.0462497a] a 
14.5032535 - 0.2644064a - 0.05037120^ - 1.5960365A^/] + 0{a^) 



= [0.6265587 + 0.1018984a] a 



+ 



7.8618848 + 1.9742471a + 0.36089030^ - 0.4799771A(/^] + 0{a^) , 



[0.7602518 + 0.1698307a] a 



+ 



8.2392595 + 2.2694139a + 0.59762970^ - 0.7417023iV/] + 0{a^) , 



— = [0.9886158 + 0.4198773a] a 



+ 



15.7278789 + 4. 1413534a + 1.5021215a^ - 1. 3862056 A(/^] + O(a^) , 



_ = [1.3116507 + 0.8520382a] a 



+ 



30.3277430 + 7.5900655a + 3.07436580^ - 2.4134871iV/] + 0{a^) , 



^f^^{p,<l)\y^a = [0.4661270 + 0.3425460a] a 



+ 



14.0528425 + 4.5189850a + 1.20622770^ - OAmmbNf 
[0.6872635 + 0.4784106a] a 

+ [20.1250083 + 5.4632742a + 1.71674350^ - 0.9665767iV/ 
[0.8137291 + 0.4321609o] a 

+ [19.0859295 + 4.8307907a + 1.5595069a2- 1.0996648A(/^]a2 + 0{a^) , 
[0.8455237 + 0.2037969a] a 

+ [l0.9356060 + 2.6215347a + 0.73451790^ - 0.8499747iV/] + 0{a^) , 
- 0.1481481a 

2.9461750 + 0.4422360a + 0.0123457a2 - 0.1476260iV/l + O(a^) , 



= - 0.3472455a 



6.6312283 + 0.0580777o + 0.02893710^ - 0.5015769A(/^] a^ + 0{a^) , 



Ef}Y.'(p,q) = - 0.5463429a 

(14) V/-' 'iJ MS 



- 10.3162815 - 0.3260806a + 0.04552860^ - 0.8555277iV/] 
+ 0{a^) 



a 



(6.3) 



and 



^(1) 



{p, q) 



MS 



-(2) 



{p,q) 



MS (3) 
- 0.3333333 + [0.5416434 - 0.1943979aJ a 

+ [2.0416107 - 1.0743003O - 0.6039997o^ + 0.0787529iV/ 

+ 0{a^) , 



MS 



^(4) 



{p, q) 



MS 



[1.4720825 + 0.3056953o] a 
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+ 



18.7974908 + 4.5957818Q + 1.0954083a^ - 1.3299518iV/ 



MS 



[1.5739809 + 0.6019916a] a 



MS 



+ 



(P, 9) 



MS 



MS 



27.3251890 + 7.1002047a + 2.1571366a2 - 1.8413671iV/ 
[1.6758793 + 0.8982879a] a 

+ [35.8528872 + 9.6046276a + 3.2188650a^ - 2.3527823A(/- 



MS 



MS 

[1.7777778 + 1.1945842a] a 

+ [44.3805855 + 12.1090504a + 4.2805934a2 - 2.8641975iV/ 

+ O(a') , 

^(13) ^^''■^^MS "(14) ^^'■^'MS 



,q) = S 

MS 

0.6944910a 



- 13.2624565 + 0.1161554a + 0.0578743a^ - 1.0031537iV/ 
+ 0{a^) . 



(6.4) 



With the two loop corrections computed we note that the results for some of the dW^ and ddWs 
amplitudes are proportional to those of lower moment operators. These relations were noted 
at one loop in [10] and extend to two loops now. For instance, channels 4 and 7 of ddWs are 
related to channels 3 and 5 respectively of dW2, [SKTO], or equally channels 2 and 3 of the vector 
current, |101lllj. Also channels 4 and 7 of dWs are proportional to channels 5 and 7 respectively 
of W2- That not all channels have a similar relation is due to the imbalance of indices of the 
operator inserted in the Green's function. Also the absence of direct equality stems from the 
difference in the two bases used for each operator moment. Though this underlying agreement 
is a partial check on our computation. Equally the equality of several of the ddW^ channels 
with themselves reflects the symmetric nature of this particular operator and is another minor 
calculational check. 

The expressions for the RI'/SMOM scheme amplitudes have a similar structure to the MS 
ones. However, given the nature of the scheme definition which we have used then there are no 
corrections to channels 1 and 2 of W3 as well as 1 of dW^. For channels 3 of and 2 and 3 of 
dWs only the tree term is present which we record with the rest of the amplitudes. We have 



^(3)' 



- 0.3333333 + O(a^) , 
[0.3425460 + 0.0484212a] a 

+ [4.8378239 + 1.2568504a + 0.1917623a2 + 0.0363159a^ 

- (0.2361911 + 0.0538013a)iV/] + O(a^) , 
[0.4238476 + 0.0751598a] a 

+ [5.1111487 + 1.4983490a + 0.2735246a^ + 0.0563697a^ 
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- (0.3660535 + 0.0835109a)Af/] o? + 0{a^) , 
[0.5842793 + 0.1503196a] a 



+ 



5.9461385 + 2.2633283a + 0.5397956a^ + 0.11273970^ 

- (0.6841872 + 0.1670218a)Ar/]a2 + O(a^) , 
1.0677459 + 0.6648679a] a 

27.4054010 + 13.2965737a + 3.36541 19a^ + 0.49865090^ 

- (2.1261604 + 0.7387421a)Arj] + 0{a^) , 
[0.2222222 + 0.1553757a] a 



+ 



+ 



6.7973293 + 3.2105607a + 0.57728610^ + 0.1165318a^ 

- (0.1633838 + 0.1726396a)iVj] + O(a^) , 
[0.2829270 + 0.2088529a] a 

+ [8.7988864 + 4.2662444a + 0.86011080^ + 0.1566397o^ 

- (0.2645584 + 0.2320588o)iV/] + O(a^) , 
[0.4773248 + 0.3374900o] a 



+ 



16.2318808 + 7.5409748O + 1.70909010^ + 0.25311750^ 

- (0.7240161 + 0.3749886a)Ar/] + 0{a^) , 
[0.5615110 + 0.1503196a] a 

+ [9.0337049 + 3.3169665a + 0.74938750^ + 0.1127397a^ 

- (0.6061888 + 0.1670218a)Arj] + 0{a^) , 

- 0.0607048a 

- 1.0472892 + 0.1315372a + 0.00505870^ - 0.0499049iV/l + 0{a 

- 0.1214095a 

- 2.2043458 + 0.1489816o + 0.01011750^ - 0.1270706A(^] a^ + 0{a 

- 0.4588995a 

- 10.4580788 - 0.5150531o + 0.03824160^ - 0.7578066iV/" 
+ 0(a3) , 



- 0.1666667 + O(a^) , T,^^\p,q) = - 0.3333333 + O(a^) 
[0.6265587 + 0.1018984a] a 



+ 



8.4968205 + 2.3224637o + 0.4090233a^ + 0.0764238a' 
- (0.4799771 +0.1 132205a) A^fja^ + 0{a^) , 
[0.7602518 + 0.1698307a] a 



+ 



9.3810330 + 2.9919683a + 0.63674050^ + 0.1273731o^ 
- (0.7417023 + 0.1887008a)A(f]a2 + 0{a^) , 
[0.9886158 + 0.4198773o] a 
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(12) 



+ [l7.8616631 + 7.6725826a + 1.8610530a^ + O.SUOOSOa^ 

- (1.3862056 + 0.4665304a)/^/]a2 + 0{a^) , 
[1.3116507 + 0.8520382a] a 

+ [33.9387109 + 16.3643066a + 4.0819609a^ + 0.6390287a^ 

- (2.4134871 + 0.9467091a)iV/] + O(a^) , 
Sfg^' (^'' 9) = [0.4661270 + 0.3425460a] a 

+ [13.3306392 + 6.3052936a + 1.2938351a2 + 0.2569095a^ 

- (0.4507105 + 0.3806067a)iV/] + O(a^) , 
^f^' {P, q) = [0.6872635 + 0.4784106a] a 

+ [20.7144109 + 9.6754987a + 2.1813682a2 + 0.3588079a^ 

- (0.9665767 + 0.5315673a)A(/^] + O(a^) , 
Sf|^'(p,g) = [0.8137291 + 0.4321609 



la a 



+ 



20.5017650 + 9.0345941a + 2.0723060a^ + 0.3241207a^ 
- (1.0996648 + 0.4801788a)A(/^] + 0{a^) , 
^{u) (^^' 9) = [0-8455237 + 0.2037969a] a 

12.6927015 + 4.3825798a + 0.9666484a2 + 0.1528477a^ 
(0.8499747 + 0.2264410a)iV/] + O(a^) , 



+ 



T,f^^p,q) = - 0.1481481a 



+ 0( 



and 



2.4900831 + 0.1788615a + 0.0123457a2 - 0.14762607Vj 

- 0.3472455a 

- 7.1954780 - 0.1444336a + 0.02893710^ - 0.5015769A(^ 

- 0.5463429a 

- 11.9008728 - 0.4677288a + 0.04552860^ - 0.8555277iV/ 
+ 0{o?) 



= - 0.3333333 + [0.5416433 - 0.2500466a] a 

+ [10.5296794 + 4.0831575a + 0.93767470^ + 0.18753490^ 

- (0.6990249 + 0.2778295o)Arj] + 0{a^) , 

= [1.4720825 + 0.3056953o] a 

+ [l8. 7974908 + 5.1040424o + 1.14635750^ + 0.2292715a^ 

- (1.3299518 + 0.3396615o)A^/]a2 + O(a^) , 
^f^^{p,q) = Sffor(p,<?) 

= [1.5739809 + 0.6019916a] a 
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+ 27.3251890 + 9.8676624a + 2.25746860^ + 0.4514937a^ 



- (1.8413671 +0.6688796a)iV/]a2 + ©(a^) , 





[1.6758793 + 0.8982879a] a 

+ [35.8528873 + 14.6312824a + 3.3685797a2 + 0.6737159a^ 



- (2.3527823 + 0.9980977a)iV/] + ©(a^) , 




[1.7777778 + 1.1945842a] a 

+ [44.3805855 + 19.3949025a + 4.4796908a2 + 0.8959382a^ 



- (2.8641975 + 1.3272318a)iV/] + ©(a^) , 






- 0.6944910a 

- 13.2624565 - 0.8098326a + 0.0578743a2 - 1.0031537iV/j 



(6.7) 



The same relations between the various channels noted earlier for the MS scheme apply to the 
corresponding amplitudes in the RI'/SMOM scheme. 

7 Discussion. 

We conclude our discussions with brief remarks. Clearly we have provided the full two loop 
structure of the Green's function with level W3 operators inserted in a quark 2-point function 
in two renormalization schemes. The underlying renormalization allowed for the computation 
of the conversion functions from the RI'/SMOM scheme to MS and we have come to the same 
conclusion as [9l |T0] that the series convergence is marginally worse than the conversion functions 
of the RI' scheme which has potential infrared issues. Though as we have argued that in some 
sense RI' is not as full a scheme as RI'/SMOM as there is no access to off-diagonal elements of the 
mixing matrix. Although one possibility of improving the convergence rests in the redefinition 
of the RI'/SMOM scheme to take account of the structure of more amplitudes. This could be 
achieved by a different choice of tensor basis for the Green's function. However, in providing 
the full structure of the Green's function at the symmetric point in MS one is free to perform 
the renormalization in any scheme of their choosing before converting to MS as the reference 
scheme. Therefore, our results are reasonably comprehensive so as to allow others to be flexible 
in how they choose to define their own version of RI'/SMOM. 



Acknowledgement. The author thanks Dr. P.E.L. Rakow for useful discussions. 

A Basis tensors and projection matrix. 

The explicit forms of the fourteen basis tensors in d-dimensions are 




{p, q) 



— [l^.PvPa + IvPaPpi + laPfji.Py\ + 



1 



[d + 2] 
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2|J 



2^ 



[d + 2]/i2 
2^ 



[d + 2]/x2 

-^PixPvPcT + _)_ 2]^2 i^M'^f''^ + ^'^o-^'m + ^o^iPv] , 



[d + 2]//2 

— IPt^QuQa + qt^PuQa + q,MquPa] 

A* 

-^PfMPuPa + +^2]/x2 ^^''''^'^ , 
[PM^'i^^o- +;>;igi/PfT + q^lPl.Pa] 

A* 

[d + 2]//^ [VfiiyPcr ~ Vfii'qa- + flvaPfi ~ Vi^a-qn "I" VcrnPi/ ~ f]<J^jiqv\ i 

A* 



1 



jqixq^qcr + _^ ^1^2 [^/'^^^^ + '^^<^^i^ + > 



+ 



r(3) npqPuPa + r(3) „pqP^xPa + ^ {^) crpgPuPv 
1 

[d + 2]/i2 

^(3) iipqPuqa + r(3) vpqP^qa + r(3) apqPf^q, 



r(3) npq^va + r(3) upqVl^a + r(3) apq'Hiii' 



1^' 



+ ^(3) ^pqqvPa + r(3) upqq^Pa + {i) apqq/^Pu 
[d + 2]/x2 -''^^^^ + ^(3) !^P?^MO- + r(3) (Tpg'?;^!/ 
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1 r. 
1 



^(3) fipgQuQa + r(3) vpqq^jLqa + r(3) apqq^Qv 



(A.l) 



We use the convention that when a Lorentz index of p^i---^" ig contracted with an external mo- 

(n) 

mentum then that momentum appears in place of the associated Lorentz index. The designation 
of in the superscript here indicates level. For dW^ and ddW'i the same basis is used. Given 
this choice we explicitly record the associated 14 x 14 matrix Ai^^, [10]. To save space we set 



M 



W3 





/ 











1 













2106d{d - 2) 




Ml' 











V 











•'^'44 



and emphasise that the r(3) sector corresponds to the lower outer corner submatrix or subspace. 
The clear partition is a reflection of the trace, (j2.18p . The explicit elements of each submatrix 
are 



'■22 



M 



23 



MY,' 



M^^ 



( 312(d + l) 156(d + l) 78(d + 4) 

156(d + l) 39(5(i + 2) lbQ{d + l] 

78(d + 4) 156(d + l) 312((i + 1 

\ 1248(d + 1) 624(d 

624(d + 1) 
312(2d + 1) 
156(3d + 4) 



1248(d + i; 
624(d + 1) 

i-ij oiZ(^a-t-ij 312(d + 4) 

+ 1) 312(d + 4) 1664(ci + 3)(d + l) y 

312((i + 2) 156(d + 4^ 

156(3(i + 2) 312(d + l 

156(3d + 4) 624(d + 1) 624(d + 1 

\ 832(d + 3)(d + l) 416(d + 6)(d+ 1) 208(d + 12)((i 

( 624((i + l) 156(3d + 4) 312((i + 

312(2d + l) 
624(d + 1 



624((i 

156(3d + 2) 

312((i + 2) 624(d+l) 1248(o! + 1) 

\ 416(d + 6)(d + 1) 208(d + 12)(d + 1) 104(^2 + 22^ + 48)/ 

/ 624(d + l) 312(2d + l) 156(3d + 4) 832(d + 3)(d + 1) \ 

312(d + 2) 156(3d + 2) 624(d + 1) 416(d + 6)(d + 1) 

156(d + 4) 312(d + l) 624((i + 1) 208(d + 12)(d + 1) 



312((i + 4) 
624((i + 1) 
1248(d + 1) 



624((i + 1) 
312((i + 1) 
156(d + 4) 
832(d + 3)(d + 



1)7 



\ 624(d + 1) 312(d + 1 

/ 416(2d + 3) 
624(d + 2) 

m/l (4d^+25d+12) 

416(d + 3) 



(d + 1 



u.i'iV"-^^; 208(d + 12)(d + 1 
156((i + 4) 832((i + 3)(d + 1) / 

624(d + 2) 







in/l (4a!2+25a!+12) 





416((i + 3) 

416(2(i + 3) 208(d + 6) 

416(4d + 3) 104(d + 12, 
104(d+12) 416(4d + 3) y 

208((i + 12)(d+ 1) \ 
416((i + 6)(d+ r 



416(2(i + 3) 
V 416(d + 3) 208(d + 6) 

/ 416(d + 3)((i + l) 312((i2 + 7d + 4) 
312((i + 2)2 416(d + 3)(d+ 1) 

208(ci + 6)((i+ 1) 416(d + 3)(d + l) 832(d + 3)((i + 1) 

\ 416(2d + 3)(d + l) 104(4^2 + 25d + 12) 208(d + 12)(d + 1) / 

/ 624(d + l) 156(3d + 2) 312(d + 2) 416(d + 6)(d + 1) ' 

156(3d + 4) 312(2d + l) 624((i + 1) 208((i + 12)((i + 1) 
\ 312(d + 4) 624(d + l) 1248(d + 1) 104(^2 + 22^ + 48)^ 

/ 416(d + 3)(d + 1) 312(d + 2)2 208(d + 6)(d + 1) 416(2d + 3)(d + 1) 
312(^2 + 7d + 4) 416(d + 3)(d + l) 416(d + 3)((i + l) 104(4(i2 + 25d + 12) 
y 208(d + 12)(d + 1) 416(d + 6)(d + l) 832(d + 3)((i + 1) 208(d + 12)(d + 1) 
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/ 208{A(f + 7d + 6) 624(d + 2) (ti + 1) 416(d + 6) (ci + 1) \ 

624(d + 2)(d+l) 416(2d + 3)(d + l) 832(d + 3)(d + 1) 
\ 416(d + 6)(d + l) 832(ci + 3)(d + l) 1664(d + 3)(d + 1) / 

/-416(d + l) -208((i + l) -104(d + 4) \ 

-208(d + l) -52(5(i + 2) -208(d + 1) . (A.2) 
^ -104(d + 4) -208(c? + l) -416(d + l) y 

These have been given in c/-dimensions because we used dimensional regularization in our two 
loop computation. 
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